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摘要 : 日 冕 物质 抛射 ( Coronal Mass Ejection, CME) 是 一 种 规模 巨大 ， 程 度 剧 烈 的 爆 
发 现象 ， 是 影响 地 球 的 主要 太阳 爆发 活动 。 由 于 这 种 爆发 现象 会 对 地 球 环境 造成 严重 干扰 ， 
因此 CME 的 探测 对 预报 灾害 性 空间 天 气 具 有 重要 意义 。 为 了 更 清楚 地 梳理 目前 存在 的 CME 


检测 方法 ， 文 章 将 对 典型 的 方法 进行 分 析 总 结 。 首 先 ， 介 绍 日 冕 物质 抛射 及 其 特征 ; 然后 ， 


从 基于 手工 方法 和 自动 检测 方法 两 方面 对 CME 检测 进行 概述 和 分 析 ; 最 后 ,讨论 目前 算法 
存在 的 一 些 问 题 ， 进 而 提出 未 来 的 研究 方向 。 

Ai: 太阳 爆发 活动 ;日冕 物质 ; CME 检测 
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70 年 代 早期 的 Skylab 任务 中 被 首次 观测 到 ， 是 源 自 太阳 的 最 大 和 最 上 其 活力 的 爆发 现象 ， 并 
且 可 以 通过 白光 日 园 仪 在 延伸 的 日 园 中 观察 到 中 。 由 于 CME 是 一 种 巨大 的 磁化 等 离子 体 喷 
发 ， 因 此 它 可 以 以 高 达 每 秒 数 千 公里 的 速度 传播 后。 鉴于 这 种 爆发 活动 引发 的 太阳 高 能 粒子 
事件 和 地 厂 暴 会 影响 航空 安全 、 卫 星 运行 、 通 信 系 统 、 电 力 设施 等 ， 因 此 有 必要 对 CME 进 
行 检测 来 为 空间 天 气 预报 人 员 提 供 预 警 信号 ， 从 而 采取 相应 的 措施 以 避免 不 必要 的 损失 中 。 
此 外 ， 有 关 CME 的 统计 信息 对 于 更 好 地 了 人 解 其 性 质 也 非常 重要 。CME 的 概念 一 直 比 较 模 
糊 ,这 种 定义 的 不 明确 性 会 导致 不 同 检测 方法 得 到 的 CME 特征 和 数量 有 所 不 同 ,.Hundhausen 


等 人 将 CME 定义 为 发 生 在 几 分 钟 到 几 个 小 时 的 时 间 尺 度 上 ， 在 日 园 视 场 中 出 现 的 新 的 ， 离 
散 的 ， 明 亮 的 白光 特征 外 。 由 于 CME 现象 与 其 他 太阳 活动 有 相似 的 特征 ， 为 了 对 CME 的 
认识 更 加 准确 ，Schwenn 等 人 将 CME 的 定义 修改 为 : 在 日 冕 视 场 中 出 现 的 一 种 新 的 ， 离 散 
的 ， 明 亮 的 白光 特征 ， 并 具有 径 向 向 外 的 速度 。 如 图 1 所 示 是 来 源 于 空间 环境 预报 中 心 
(NSSC) 中 的 SOHO/LASCO 卫星 观测 到 的 一 次 日 冕 物质 抛射 的 爆发 过 程 。 
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图 1 SOHO LASCO C2 观测 到 的 日 冕 物质 抛射 的 爆发 过 程 
Fig. 1 ACME observed by SOHO LASCO C2 

随 着 SOHO, Wind 和 STEREO 等 卫星 的 发 射 , CME 的 检测 和 识别 成 为 可 能 。 其 中 1995 
年 发 射 的 SOHO 卫星 搭载 的 LASCO 日 冕 仪 可 以 观测 到 CME 的 发 生 过 程 ， 得 到 CME 的 影 
像 资料 ， 通 过 这 些 资料 人 们 可 以 检测 CME 的 运动 状态 ， 进 而 得 到 CME 的 描述 参数 〈 位 置 
角 、 角 宽度 、 速 度 等 )。 目 前 国内 外 关于 CME 的 检测 研究 已 经 取得 很 大 进展 ， 形 成 了 许多 
CME 目录 ， 为 以 后 的 进一步 研究 提供 了 基准 和 参考 。 

为 了 方便 读者 详细 了 解 目 前 CME 检测 的 研究 进展 ， 并 为 相应 研究 提供 借鉴 ， 本 论文 梳 
理 了 国内 外 关于 CME 检测 的 一 些 方法 ， 给 出 了 具有 代表 性 方法 的 具体 识别 过 程 。 同 时 为 了 
便于 读者 理解 ， 本 论文 将 目前 CME 检测 方法 分 成 了 两 大 类 : 基于 手工 识别 的 方法 和 基于 自 
动 检测 的 方法 , 同时 又 根据 所 选取 的 特征 和 采用 的 技术 将 自动 检测 方法 分 为 四 类 : 基于 灰 度 
特征 、 基 于 纹理 特征 、 基 于 光 流 法 和 基于 学 习 的 方法 来 详细 概述 。 

本 文 总 共 分 为 四 部 分 : 第 一 部 分 为 引言 ， 引 出 CME 的 定义 以 及 CME 检测 的 重要 意义 ; 
第 二 部 分 介绍 CME 的 特征 ， 包 括 速 度 、 张 角 和 发 生 规律 ， 第 三 部 分 分 类 总 结 目前 现 有 的 
CME 检测 方法 ， 第 四 部 分 总 结 与 展望 ， 总 结 目 前 CME 检测 存在 的 挑战 ， 探 讨 未 来 的 研究 
趋势 。 本 文 的 贡献 在 于 ; 

C1) 目前 未 见 有 相似 文献 ， 可 为 研究 工作 者 提供 该 领域 的 研究 现状 ; 

(2) 可 供 研 究 工作 者 详细 了 解 利 用 不 同方 法 进行 CME 检测 的 过 程 ; 

(3) 针对 目前 CME 检测 存在 的 挑战 ， 提 出 未 来 的 研究 趋势 ， 为 研究 工作 者 的 研究 方 
向 起 到 借鉴 作用 。 


2 CME 特征 


CME 检测 实质 上 就 是 利用 人 工 标记 或 计算 机 自动 识别 的 方法 检测 CME 事件 从 而 得 到 
CME 特征 。 常 用 的 CME 特征 包括 : 速度 、 张 角 、 发 生 规律 等 。 
2.4. CME 的 速度 

CME 是 一 种 运动 的 过 程 ， 其 速度 是 一 个 重要 的 物理 特性 。 由 于 目前 观测 到 的 CME 通 
常 是 其 在 天 空 平面 的 投影 ， 因 此 一 般 情 况 下 CME 的 速度 是 指 CME 最 快 前 沿 的 运动 速度 中 ， 
它 是 由 CME 前 沿 的 “时 间 - 高 度数 据 经 过 线性 拟 合 得 到 的 。 如 图 2 是 来 源 于 协调 数据 分 析 研 
讨 会 中 心 (CDAW') 目录 中 2017 年 7 月 2 日 13 时 25 分 42 秒 发 生 的 一 个 CME 的 时 间 - 高 
FER. CME 的 速度 范围 为 10-3500 km /s, 一 般 速 度 >1000kmys 的 为 快速 CME, 速度 <100km/s 
的 为 缓慢 CME, 
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2 CDAW 目录 中 2017 年 7 月 2 号 13 时 25 分 4 和 2 秒 发 生 的 一 个 CME 的 时 间 一 高 度 图 


Fig. 2 The time-height plot of a CME occurring at 13:25:42 on July 2, 2017 in the CDAW catalog 


22 CME 的 张 角 

CME 的 大 小 一 般 用 张 角 《〈 角 宽度 ) 来 衡量 ， 张 角 是 CME 区 域 两 个 边缘 所 构成 的 角度 ， 
在 数值 上 等 于 两 边缘 的 位 置 角 之 差 。 图 3 是 来 源 于 CACTus 软件 检测 到 的 一 个 CME， 其 中 
a. b 两 处 为 某 时 刻 发 生 的 CME 区 域 的 位 置 角 ， 张 角 则 为 图 中 白色 实 线 中 间 的 角度 。CME 
的 张 角 从 几 度 到 360 度 ，360 度 的 CME 称 为 Halo CME， 一 般 宽度 >120 度 的 为 宽 CME, Sit 
度 <30 度 的 为 窗 CME", 


2000/01/02 07:54:05 


图 3 CME 的 角 宽 度 


Fig. 3 The angular width of a CME 


2.3 CME 与 太阳 黑子 数 的 关系 

CME 的 发 生 规律 与 太阳 黑子 活动 周 有 关 ， 在 太阳 活动 极 小 年 时 为 0.5 个 /天 ， 在 太阳 活 
动 极 大 年 时 为 2-6 个 /天 。 下 图 4 表示 2008 年 到 2017 年 间 CDAW 目录 记录 的 CME 数量 ， 
5 中 的 数据 是 来 源 于 太阳 黑子 指数 和 长 期 太阳 观测 CSILSO?O. 中 的 近 13 年 太阳 黑子 的 数 
量 , 我 们 可 以 看 到 在 这 两 个 图 的 重茬 时 间 段 (2008 年 一 2017 年 ) CME 的 发 生 频 率 和 黑子 的 
基本 一 致 ， 因 此 可 以 说 明 CME 的 发 生 规 律 与 太阳 黑子 活动 周 有 关 。 
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图 4 2008 年 到 2017 年 间 CDAW 目录 记录 的 CME 数量 
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International sunspot number S, : last 13 years and forecasts 
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图 5 近 13 年 来 太阳 黑子 数量 变化 情况 


Fig. 5 Evolution of sunspot numbers in the past 13 years 


3 CME 检测 方法 


十 多 年 来 ， 研 究 者 们 已 经 利用 多 种 方法 检测 到 CME 事件 及 其 特征 (例如 速度 ， 角 宽度 ， 
发 生 规律 等 )， 并 且 制 成 了 CME 目录 来 供 科研 人 员 学 习 研 究 。 目 前 CME 的 检测 方法 主要 分 
为 两 大 类 : 

(1) 基于 手工 方法 标识 CME : 目前 主要 有 两 个 著名 的 目录 : CDAW 目录 和 海军 研究 实验 
室 (NRL) 目录 , 这 种 方法 主要 是 靠 观 察 者 每 天 手动 记录 数据 来 进行 CME 编目 ， 比 较 费 时 ; 
(2) 基于 自动 方法 检测 CME: HF CME 是 一 种 动态 变化 的 过 程 ， 通 常 具有 明亮 的 强度 特 
征 和 复杂 的 纹理 结构 。 我 们 可 以 依据 CME 检测 选用 的 特征 及 技术 手段 将 目前 现 有 的 自动 检 
测 方法 分 为 以 下 四 类 来 进行 概述 和 分 析 ， 如 下 表 1 所 示 ， 是 对 目前 CME 检测 方法 的 总 结 。 
表 1 CME 检测 方法 总 结 


Tab.1 Summary of CME detection methods 


CME 检测 方法 


NRL 


CDAW 


Berghmans ^: A. (2002, 2004) P?! 
Olmedo 等 人 (2008, 2010) ° 
Qu ^£ (2006) | 


Boursier “£ A (2005, 2009) U? 5! 


基于 灰 度 特征 方法 


自动 检测 方法 


Gonzalez-Gomez 等 人 (2010) 0 
Gallagher 等 人 (2011) 5 

基于 纹理 特征 方法 Goussies 等 人 (2010) 05 171 
Braga 等 人 (2013) H% 


Morgan 等 人 (2012) 091 
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曾 昭 宪 等 人 (2012) P 


Colaninno 等 人 (2006) P" 
基于 光 流 法 
Gissot 等 人 (2003) °” 


Zhang 4& A (2016, 2017) 74! 


基于 学 习 方 法 Yin 等 人 (2017) P?l 
姚 海 等 人 (2017) P! 
3.1 手工 识别 目录 与 自动 检测 目录 的 比较 
CME 与 耀 班 ， 太 阳 高 能 粒子 ， 地 磁 风 暴 等 许多 现象 有 关 ， 因 此 检测 CME 现象 并 编制 
事件 目录 对 于 国内 外 研究 工作 者 理解 这 些 现象 并 开展 相关 工作 非常 重要 。CME 目录 是 将 检 
测 到 的 CME 事件 特征 以 数据 和 图 表 的 形式 记录 下 来 ， 相 关 参 数 特征 主要 包括 : 速度 、 角 宽 
度 、 位 置 角 等 。 根 据 上 述 介绍 的 CME 检测 方法 我 们 可 以 得 到 两 类 CME 目录 : 手工 识别 目 
录 和 自动 检测 目录 ， 通 过 以 CDAW 目录 和 CACTus 目录 为 代表 进行 特征 〈 速 度 和 角 宽 度 ) 
的 比较 四， 我 们 得 到 了 这 两 类 目录 的 特点 : 
C1) 识别 狭窄 的 CME CW«300: 研究 表明 ， 在 同一 时 期 CACTus 检测 到 的 罕 CME 远 
远 多 于 CDAW， 主 要 原因 在 于 CDAW 目录 在 太阳 活动 极 大 年 遗漏 了 许多 窄 CME. 
(2) 识别 快速 的 CME (V>1000km/s): 在 同一 时 期 虽然 CACTus 检测 到 的 快速 CME 
多 于 CDAW， 但 是 在 检测 到 的 32% 的 真正 CME 中 只 有 6% 是 快速 CME. 
因此 可 以 得 出 结论 : 手工 目录 和 自动 检测 目录 各 有 优 人 缺点， 由 于 手工 目录 是 基于 人 标记 
的 ， 因 此 会 不 可 避免 的 遗漏 许多 狭窄 的 CME; 而 自动 检测 目录 虽然 比 手工 目录 识别 的 快速 
CME 多 ， 但 是 准确 率 不 太 好 。 
3.2 手工 方法 标识 CME 
自从 1971 年 发 现 CME 现象 以 来 , CME 已 经 被 多 种 星 载 仪 器 观测 到 。 自 1995 年 SOHO 
航天 器 发 射 以 来 ，CME 观测 主要 采用 LASCO HZA. LASCO CME 的 识别 和 编目 是 一 项 
重要 任务 ， 为 进一步 的 科学 研究 提供 基础 知识 。 
NRL 目录 由 LASCO 观察 员 编制 , 他 们 查看 LASCO 日 铣 图 像 的 序列 并 记录 每 天 发 生 的 
事件 。 这 是 一 份 初步 目录 ， 提 供 有 关 CME 时 间 和 大 致 位 置 的 信息 。 
CDAW 目录 是 由 美国 宇航 局 、 美 国 天 主教 大 学 以 及 海军 研究 实验 室 合作 ， 在 CDAW 数 
据 中 心 生成 和 维护 。 标 记 的 每 个 CME 由 发 生 的 日 期 和 时 刻 确定 ， 该 目录 提供 了 CME 属性 
的 一 些 测量 值 ， 包 括 速 度 和 角 宽 度 等 ， 同 时 还 提供 了 确定 CME 速度 的 “时 间 - 高 度 ” 图 。 除 了 
© 得 到 的 数据 和 图 表 之 外 , 该 目录 还 包含 由 日 冕 仪 得 到 的 CME 图 像 序列 , 可 供用 户 下 载 观 看 。 
由 于 NRL 目录 和 CDAW 目录 都 是 由 观察 者 手工 标识 的 , 即使 对 于 同一 观察 者 , 识别 能 
力也 不 是 恒定 的 。 因 此 ,所 提供 的 事件 和 测量 的 参数 受 人 类 主观 因素 影响 较 大 ， 同 时 编目 过 
程 也 非常 耗 时 。 
3.3 自动 方法 检测 CME 
传统 的 CME 检测 方法 是 基于 人 类 观察 ， 这 种 观察 方法 效率 低 ， 容 易 受 到 个 体 主观 因素 
的 影响 ， 而 随 着 近年 来 自动 识别 技术 的 迅速 发 展 ， 出 现 了 很 多 自动 检测 CME 的 算法 ， 为 
CME 的 检测 研究 提供 了 新 方法 。 由 于 CME 具有 不 同 的 形态 ， 尺 度 和 特征 ， 通 常 表现 为 一 
个 明亮 的 、 纹 理 复杂 的 增强 结构 ， 尾 随 着 一 个 亮度 不 足 的 暗 区 域 申 ， 我 们 可 以 将 目前 现 有 的 
自动 CME 检测 方法 分 为 以 下 四 类 。 
3.3.1 基于 灰 度 特 征 的 方法 
灰 度 特征 是 基于 像素 点 表示 图 像 亮度 的 特征 ， 与 背景 灰 度 特征 不 同 ，CME 具有 更 明亮 
的 结构 ， 通 常 表现 出 明亮 的 白光 特征 ， 因 此 通过 提取 灰 度 特征 检测 CME 是 一 种 很 直接 的 思 
路 。 目前， 基于 灰 度 特征 检测 CME 的 方法 研究 已 较 成 熟 ， 国 内 外 已 经 有 不 少 人 利用 该 方法 
进行 CME 检测 并 且 取 得 了 很 好 的 效果 。 计 算 机 辅助 跟踪 软件 包 CCACTus) 和 太阳 能 爆发 
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事件 检测 系统 CSEEDSO 是 其 中 具有 代表 性 的 两 种 检测 方法 。 

CACTus 方 法 由 Berghmans 等 人 W' | 提 出 ， 该 方法 通过 滤波 去 品 、 极 坐标 变换 、 帧 间 差分 
以 及 整合 LASCO C2 和 C3 图 像 完 成 预 处 理 。 对 不 同时 刻 t 的 差分 图 像 [9,7] 进 行 整 合 ， 得 到 [z， 
0, 7] 数据 立方 体 , 该 立方 体 沿 指 定 角 度 9 的 切片 [t, J] 中 的 一 条 次 线 是 一 个 CME 沿 指定 角度 9 
上 的 喷射 表现 。 为 了 能 清楚 的 看 到 疹 线 ,将 切片 图 颜色 进行 反 转 得 到 倾斜 的 暗黑 疹 线 ， 如 图 
6 所 示 。 由 于 霍 夫 变 换 C9 是 一 种 图 像 处 理 技术 ， 它 可 以 从 图 像 中 检测 出 线段 ， 因 此 霍 夫 变 换 
被 用 来 检测 [5， 门 切片 中 的 兰 线 。 同 时 霍 夫 变换 检测 到 的 兰 线 是 沿 一 个 角度 切片 的 CME， 为 
了 得 到 完整 的 CME， 将 图 像 投 影 到 [v，69， 避 数据 立方 体 中 ， 速 度 v 可 通过 兰 线 段 长 度 计算 
得 到 ， 此 时 CME 识 别 便 转换 为 [vy，0， 可 立方 体 中 点 的 聚 类 问题 。 由 于 CACTus 认 为 同一 个 
CME， 治 不 同 g 方 向 的 径 向 速度 相似 ，v 方 向 信息 对 聚 类 贡献 较 小 ， 为 了 降低 计算 复杂 度 ， 
对 [v，0, 本 沿 vy 方 向 积分 ， 使 得 聚 类 转换 到 二 维 空间 [9， 中 进行 从 而 在 CME 概 览 图 中 得 到 
聚 类 的 位 置 ， 该 位 置 为 CME 发 生 的 起 始 时 间 和 和 角度， 垂直 方向 上 聚 类 的 长 度 表 示 CME 喷 发 
持续 的 时 间 。 
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图 6 时 间 - 高 度 图 中 检测 次 线 ， 图 像 来 源 于 Berghmans 等 (2004) 加 


Fig. 6 Detected ridges in a time-height diagram, Image from Berghmans et al. (2004) ?! 

CACTus 软件 是 第 一 个 利用 自动 方法 实现 CME 检测 ， 相 比 于 人 工 标识 ， 该 方法 检测 速 
度 快 ， 检 测 准确 率 也 很 好 ， 可 以 检测 到 CDAW 目录 中 约 75% 的 CME， 并 且 也 能 够 检测 到 
CDAW 目录 中 没有 检测 到 的 弱 CME。 目前 , CACTus 软件 编制 了 一 个 从 1996 年 到 现在 的 在 
HC. 

SEEDS 方法 由 Olmedo 5$ AP 四 提出 ,检测 过 程 使 用 基于 区 域 增长 算法 的 图 像 分 割 技术 。 
与 CACTus 类 似 ， 该 方法 通过 归 一 化 输入 图 像 、 滤 波 去 品 、 极 坐标 变换 完成 预 处 理 。 然 后 将 
极 坐 标 变 换 得 到 的 [9,7] 图 像 序列 沿 0 方向 进行 投影 得 到 信和 号 强度 曲线 , 如 图 7(a) 所 示 。 接着 
通过 闵 值 处 理 得 到 CME 的 中 心 和 最 亮 的 部 分 ， 即 核心 角度 ， 同 时 对 核心 角度 利用 区 域 增 长 
算法 得 到 整个 CME 的 角 宽 度 。 CME 前 沿 的 确定 是 通过 将 图 像 沿 7 方 向 进行 投影 得 到 信号 强 
度 曲线 ， 如 图 7(b) 所 示 ， 曲 线 峰值 对 应 于 Max-Height 处 ， 曲 线 峰值 的 一 半 在 Half-Max-Lead 

CATH) U Half-Max-Follow〔 后 继 ) 处 。 通 过 对 比 连续 两 帧 图 像 ， 如 果 后 张 图 像 的 前 沿 和 
后 继 以 及 两 者 之 差 高 于 前 张 图 像 时 ， 表 示 CME 是 扩张 的 ， 同 时 车 前 沿 溢出 视 场 或 者 CME 
太 瞳 识别 不 出 时 ， 就 会 停止 追踪 ， 寻 找 下 一 个 CME。 


3 CACTus: http://sidc.oma.be/cactus/ 
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图 7(a) 沿 6 轴 的 一 维 投影 强度 分 布 图 ， 图 7(b) 沿 r 轴 的 一 维 投影 强度 分 布 图 ， 
图 像 来 源 于 Olmedo 5 (2008) P! 图 像 来 源 于 Olmedo & (2008) P! 
Fig. 7(a) The intensity profile of the 1D projection along Fig. 7(b) The intensity profile of the 1D projection along 
the 0 axis, Image from Olmedo et al. (2008) 中 the r axis, Image from Olmedo et al. (2008) 四 


在 检测 准确 率 上 ， 与 CACTus 类 似 ，SEEDS 能 检测 到 CDAW 目录 约 7590] CME. {A 
在 检测 到 的 CME 总 数量 上 ， 通 过 CACTus 和 SEEDS 方法 检测 得 到 的 CME 数量 是 CDAW 
目录 的 两 倍 多 。 同 时 SEEDS 方法 还 可 以 得 到 CME 的 运动 前 沿 轮廓 ， 这 可 以 在 SEEDS 目录 
“看 到 。 

除了 上 述 详细 介绍 的 两 种 方法 之 外 ， 还 有 许多 基于 灰 度 特征 检测 CME 的 方法 。 例 如 
Qu 等 人 0 提出 了 一 种 基于 支持 向 量 机 〈Support Vector Machine, SVM) 的 CME 检测 方法 ， 
为 了 将 CME 与 其 他 结构 区 分 开 来 ， 该 算法 通过 预 处 理 产 生 差分 图 像 和 分 割 图 像 ， 然 后 再 次 

利用 阔 值 方法 来 分 割 预 处 理 得 到 的 差分 图 像 和 分 割 图 像 , 最 终 的 分 割 结果 是 来 自 两 个 分 割 结 

果 的 总 和 。 接着 将 所 有 分 割 区 域 视 为 CME 候选 区 域 , 通过 设 定 阔 值 将 候选 区 域 分 类 为 CME 
区 域 和 背景 ， 最 后 在 CME 区 域 利用 SVM 分 类 器 来 区 分 强 CME 4159 CME。 与 其 他 方法 相 
比 ， 该 方法 的 最 大 优点 是 能 够 对 强 CME 和 大 多 数 弱 CME 进行 准确 而 快速 的 检测 。 

与 上 述 三 种 方法 将 LASCO 图 像 转换 为 差分 处 理 的 极 坐 标 图 不 同 ，Boursier 等 人 提出 的 
ARTEMIS FYR” VY LASCO C2 图 转换 为 天 气 图 。 如 图 8 所 示 ， 天 气 图 的 横 坐 标 为 时 间 ， 
表示 CME 的 流出 时 间 ， 纵 坐标 为 纬度 ， 表 示 CME 的 角度 宽度 。 该 算法 利用 中 值 滤波 器 去 
除 噪声 ， 然 后 对 天 气 图 进行 阔 值 处 理 返回 一 个 二 维 撼 码 ， 该 掩 码 定义 了 图 像 的 感 兴趣 区 域 。 
最 后 通过 引入 CME 特性 的 先 验 知识 正确 识别 CME。 算法 最 终 得 到 一 个 ARTEMIS 目录 , 列 
出 了 每 个 检测 到 的 CME 事件 以 及 它们 的 主要 参数 。 与 CDAW 目录 进行 比较 可 以 发 现 ， 该 
方法 的 优点 是 能 够 检测 到 一 些小 的 , 微弱 的 CME, 与 CACTus、SEEDS 一 样 , 检测 到 的 CME 
事件 多 于 CDAW 目录 。 


4 SEEDS: http://spaceweather.gmu.edu/seeds/ 
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Time 
图 8 将 CME 图 转换 为 天 气 图 ， 图 像 来 源 于 Boursier 等 (2009) M! 
Fig. 8 Converting a CME image to a synoptic map, Image from Boursier et al. (2009) !?! 


基于 灰 度 特征 的 四 种 方法 先 对 LASCO 图 像 预 处 理 得 到 极 坐标 图 或 天 气 图 ， 再 进行 阔 值 


处 理 和 形态 学 处 理 提 取 CME 特 征 ， 最 后 通过 对 CME 特 征 的 识别 或 跟踪 计算 出 CME 的 各 种 属 
TE. 我 们 可 以 发 现 这 些 方法 有 一 个 共同 点 , 那 就 是 都 采用 了 亮度 增强 的 方式 突出 日 蝎 图 像 中 
的 感 兴趣 区 域 即 CME 区 域 。 同 时 上 述 几 种 方法 在 同一 时 期 的 LASCO 图 像 上 检测 出 的 CME 数 
量 都 比 CDAW 目 录 多 , 但 也 都 不 能 完全 检测 出 CDAW 中 100% 的 CME, 这 与 自动 检测 方法 使 用 
的 图 像 预 处 理 技术 和 检测 规则 、 阔 值 的 选取 以 及 CME 的 不 明确 定义 都 有 关系 。 
3.3.2 基于 纹理 特征 的 方法 

CME 不 仅 具 有 明显 的 亮度 结构 ， 而 且 还 有 复杂 的 纹理 结构 。 对 于 一 些 比较 暗 但 具有 明 
显 纹理 结构 的 CME 区 域 来 说 ， 利 用 纹理 特征 检测 是 一 个 不 错 的 选择 。 纹 理 特征 的 提取 一 般 
需要 设 定 一 定 大 小 的 , 包含 多 个 像素 点 的 窗口 区 域 , 然后 从 中 取得 纹理 特征 ， 尤 其 在 检索 具 
有 明显 粗细 ， 玻 密 等 图 像 时 有 恨 好 的 表现 。 

Gonzalez-Gomez 等 人 使 用 小 波 分 析 方 法 对 CME 进行 分 类 , 该 算法 是 一 种 基于 频 域 分 
析 的 方法 。 通 过 在 LASCO C2 图 像 上 放置 固定 像素 的 窗口 (实验 数据 表明 ， 窗 口 可 以 放置 
在 图 像 的 任何 位 置 )， 使 用 标准 快速 傅 里 时 变换 算法 在 窗口 上 执行 卷 积 ， 同 时 使 用 墨西哥 由 
作为 基 函 数 来 获得 窗口 区 域 的 小 波谱 ， 得 到 的 小 波谱 图 横 坐 标 为 尺度 ， 纵 坐标 为 通 量 。 窗 口 
区 域 表 示 的 曲线 都 有 一 个 拐点 , 该 点 区 分 了 通 量 随 尺 度 变 化 的 速度 , 通过 计算 每 条 曲线 的 二 
阶 导 数 获得 拐点 的 尺度 (a.〉 和 通 量 (f.)。 如 图 9 R, E acf, 图 中 可 以 看 到 小 波谱 明显 被 
分 为 两 组 : 低 通 量 、 小 空间 尺度 组 (Homogeneous Group) 和 高 通 量 、 大 空间 尺度 组 (Collimated 
Group)， 其 中 box1 和 box2 表示 图 像 上 的 两 个 窗口 。 将 该 方法 应 用 到 更 大 的 CME 图 像样 本 
中 ， 通 过 在 ae- 刀 图 中 绘制 抛 点 的 尺度 和 通 量 值 ， 将 这 些 点 进行 线性 最 小 二 乘 拟 合 ， 就 可 以 
分 类 出 CME 图 像 属于 哪 一 组 。 该 方法 是 一 种 图 像 滤波 的 分 析 方 法 ， 利 用 频 域 的 高 通 和 低 通 
特征 来 分 类 CME 图 像 ， 可 以 对 由 不 同 卫 星 获 得 的 日 冕 物质 抛射 进行 直接 分 类 。 
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图 9 CME 图 像 分 类 结果 ， 图 像 来 源 于 Gonzalez-Gomez “ (2010) "^! 


Fig. 9 Classification results of CMEs images, Image from Gonzalez-Gomez et al. (2010) "4! 

CME 图 像 特征 有 多 尺度 特性 ， 多 尺度 图 像 处 理 技 术 对 增强 CME 前 治 的 可 见 性 和 抑制 
噪声 有 很 大 作用 。 由 于 线段 更 能 反应 出 图 像 信 息 ， 而 小 波 更 适合 识别 点 状 特征 ， 如 噪声 或 背 
景 恒星 ， 不 适合 检测 CME 特征 的 线 结构 ， 为 此 Gallagher 等 人 5 研究 了 高 阶 多 尺度 技术 ， 
如 痊 波 和 曲线 波 。 与 小 波 变 换 不 同 , 疹 波 变换 首先 进行 Randon 变换 ， 用 线 参 数 取 代 点 参数 ， 
之 后 进行 小 波 变 换 。 该 变换 能 够 有 效 的 表示 直线 奇异 性 特征 , 但 是 不 能 很 好 的 表示 图 像 边 缘 
曲线 。 而 曲线 波 变换 弥补 了 将 波 变换 的 不 足 , 该 方法 能 够 最 优 的 表示 图 像 的 曲线 奇异 性 特征 。 
小 波 、 关 波 和 曲线 波 类 似 ， 都 采用 基 函 数 与 信号 的 内 积 实现 信号 的 稀 玻 表示 ,不同 的 是 兰 波 
和 曲线 波 变 换 的 去 噪 效 果 好 , 同时 能 够 较 好 的 表达 图 像 边缘 信息 。 该 论文 分 别 利用 小 波 变 换 
和 曲线 波 变 换 过 滤 原 始 图 像 ， 对比 实验 结果 表明 , 与 小 波 相 比 ， 曲 线 波 变换 能 够 获得 较 好 的 
去 噪 效 果 ， 有 效 地 增强 了 CME 图 像 的 前 沿 结构 。 

图 像 纹理 可 以 表示 为 像素 灰 度 值 的 空间 变化 的 函数 , 纹理 特征 的 提取 除了 利用 信和 号 处 理 
领域 的 线性 多 尺度 变换 方法 增强 图 像 的 边缘 信息 来 检测 CME 前 治 结构 之 外 ， 另 一 种 有 效 的 
方法 是 利用 砍 度 共生 甜 阵 " 捕 获 并 表征 不 同 区 域 的 纹理 信息 。 Goussies 等 人 提出 了 一 种 基于 
灰 度 共生 和 矩阵 的 非 参 数 监督 的 CME 分 割 方法 I* 77, gy CORSET 算法 (CORonal 
Segmentation Technique)， 该 方法 是 在 有 监督 的 区 域 竞 争 模型 分 割 方法 上 的 改进 。 区 域 竞 争 
的 目标 是 将 图 像 分 割 成 多 个 区 域 , 使 得 每 个 区 域 中 的 点 具有 相似 的 图 像 特征 。 而 将 图 像 每 个 
像素 的 灰 度 值 用 作 区 域 竞 争 方法 中 的 特征 向 量 时 ， 有 监督 的 区 域 竞争 模型 分 割 CME 会 存在 
两 个 问题 : 一 是 CME 区 域 和 背景 的 灰 度 直方 图 有 重合 ， 二 是 CME 事件 和 背景 的 灰 度 直方 
图 不 遵循 正 态 分 布 ， 即 统计 模型 未 知 ， 这 两 个 问题 将 会 导致 CME 分 类 错误 。 针 对 问题 一 ， 
为 了 解决 直方 图 重 登 ， 新 算法 利用 灰 度 共生 矩阵 描述 CME 纹理 信息 ; 针对 问题 二 ， 使 用 非 
参数 模型 中 的 卡 方 统计 检验 , 它 可 以 评估 观察 到 的 事件 是 否 遵循 菜 个 特定 分 布 。 通 过 修改 区 
域 竞 争 运动 方程 ， 将 卡 方 检 验 和 纹理 信息 引入 方程 中 ， 同 时 在 图 像 分 割 中 ， 使 用 快速 水 平 集 
算法 来 实现 分 割 曲 线 的 演变 ,该 分 割 算法 首先 围绕 C2 遮挡 盘 形 成 环 状 区 域 ,如 图 10(a) 所 示 ， 
如 果 存 在 CME， 曲 线 演变 为 CME 轮廓 ， 如 图 10(b-d) 所 示 ; 如 果 不 存 在 CME， 轮 廓 就 会 消 
失 如 图 11 所 示 。 为 了 跟踪 CME, 使 用 当前 图 像 的 分 割 结果 作为 下 一 图 像 的 初始 轮廓， 并 根 
据 修 改 的 区 域 竞 争 运动 方程 演变 轮廓 ， 如 果 轮 廓 消失 意味 着 某 个 CME 事件 结束 。 该 算法 使 
用 卡 方 检验 和 灰 度 共生 矩阵， 以 便 正确 捕获 在 差分 图 像 中 观察 到 的 CME 纹理 信息 ， 可 以 检 
测 和 跟踪 具有 不 同形 状 和 强度 的 CME, 同时 得 到 了 一 个 类 似 CACTus 的 CME 概览 图 。 


a 


图 10 CME 出 现时 的 轮廓 演变 ， 图 11 没有 CME 时 的 轮廓 演变 ， 
图 像 来 源 于 Goussies 等 (2010) "7! 图 像 来 源 于 Goussies 等 (2010) "7! 
Fig. 10 Contour evolution with a CME , Fig. 11 Contour evolution without CMEs, 
Image from Goussies et al. (2010) "7! Image from Goussies et al. (2010) l” 


虽然 CORSET 算法 提供 了 CME 事件 的 边界 信息 ,但 是 没有 对 CME 参数 进行 明确 的 定 
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量 评 估 ， 为 此 在 Goussies 等 人 的 工作 基础 上 ，Braga 等 人 1 通过 添加 几 个 新 的 功能 ， 即 自动 
计算 不 同 的 形态 学 和 运动 学 参数 扩展 CORSET， 增 强 了 算法 的 功能 ， 并 将 获得 的 参数 中 
心 位 置 角 ， 角 宽度 和 速度 ) 与 现 有 的 手工 目录 和 自动 检测 目录 进行 比较 。 

由 于 CME 是 一 种 不 断 运 动 的 状态 ， 因 此 可 以 通过 将 CME 结构 与 背景 区 分 开 来 进行 
CME 检测 ， 目 前 大 多 数 研究 采用 基于 差分 的 方法 来 检测 图 像 的 运动 区 域 ， 但 是 这 种 数值 差 
分 会 将 噪声 增强 到 与 目标 相当 的 水 平 ， 虽 然 利 用 中 值 滤波 器 会 抑制 噪声 ， 但 也 会 平滑 小 型 
CME 特征 ， 同 时 差分 方法 还 会 引入 时 空 串扰 的 问题 。 为 此 Morgan 等 人 WM 提出 了 一 种 反 卷 
积 方法 ， 该 算法 利用 归 一 化 径 向 渐变 滤波 器 (NRGF) “3 将 CME 图 像 分 离 成 静态 的 背景 结 
构 和 动态 的 CME 运动 结构 ， 实 验 结果 表明 ， 该 方法 能 够 很 好 的 检测 出 弱 CME. 

CME 区 域 与 背景 以 及 其 他 太阳 结构 如 园 流 等 在 形态 、 纹 理 等 特征 上 有 明显 差异 ， 在 视 
WE, SUA 分 布 相对 均匀 ， 而 CME KRIE 22 "E BRE Jy ES ESE AO 
提出 了 一 种 基于 频谱 突变 分 析 的 CME 识别 方法 。 该 算法 与 上 述 多 尺度 变换 方法 一 样 ， 也 是 

种 基于 频 域 分 析 的 方法 , 算法 利用 傅 里 时 变换 得 到 预 处 理 后 的 日 铣 图 像 的 频谱 图 ,分离 出 
非 突 变 和 突变 信息 后 ， 通 过 傅 里 叶 反 变换 得 到 突变 结构 在 原 图 中 的 对 应 信息 来 确定 CME 的 
初步 分 离 ， 最 后 利用 局 部 稳定 极 值 区 域 检测 方法 确定 CME 区 域 的 轮廓 。 由 于 CME 区 域 轮 


= 


pu 亡 不 规则 ， 因 此 提出 基于 区 域 协 方差 的 外 接 椭圆 方法 得 到 角 宽 度 、 速 度 等 参数 。 由 于 该 方法 
o 4) TES E E ER RUE, Dott VU] d EM CME 时 具有 良好 的 效果 。 

= 上 述 CME 检测 方法 是 建立 在 时 域 、 频 域 以 及 多 尺度 分 析 的 基础 上 ， 通 过 信和 号 处 理 中 的 
c 小 波 、 兰 波 和 曲 波 等 线性 多 尺度 变换 方法 、 研 究 图 像 的 频谱 信息 以 及 利用 统计 中 的 灰 度 共生 
e 矩阵 描述 纹理 特征 ， 该 关 方 法 尤其 在 检测 弱 CME 以 及 蜡 CME 上 表现 良好 。 

©) 3.3.3 基于 光 流 法 

Ta 光 流 是 测量 数字 图 像 运动 的 一 种 强大 的 图 像 处 理工 具 ， 它 包含 图 像 中 目标 运动 的 信息 。 


光 流 算法 可 以 从 连续 图 像 序列 中 估计 出 每 个 像素 的 速度 矢量 ， 从 而 形 成 该 图 像 的 运动 场 , 可 
被 用 来 确定 目标 的 运动 情况 。 

Colaninno 等 人 提出 了 一 种 基于 光 流 法 的 CME 检测 和 跟踪 算法 ， 对 交流 约束 方程 
(OFCE) 添加 平滑 正则 化 项 马尔 科 夫 随机 场 。 然 后 将 光 流 的 估计 表示 为 全 局 优化 问题 ， 为 
了 找到 全 局 最 小 值 , 采用 多 重 网 格 松 弛 法 9 进行 求解 。 将 该 光 流 估 计 方 法 应 用 于 LASCO C2 
图 像 ， 观 察 到 1999 年 到 2004 年 间 大 约 330 个 CME 并 对 其 进行 编目 ， 最 终 对 10 个 具有 


= 显 结构 的 CME 进行 分 析 ， 推 导出 它们 的 速度 场 。 结 果 如 图 12 所 示 ， 上 面 一 行 表示 输入 


R) Zn 


Oo 像 ， 下 面 一 行 表示 光 流 结果 ， 阴 影 颜 色 表示 速度 大 小 ， 箭 头 方向 表示 在 8x8 网 格 上 采样 的 
速度 方向 。 得 到 的 速度 测量 结果 能 够 可 视 化 CME 等 离子 体 的 演化 。 该 算法 速度 非常 快 ， 可 
以 很 容易 地 应 用 于 所 有 可 用 的 CME 图 像 ， 同 时 基于 光 流 算法 检测 到 的 图 像 可 以 很 容易 地 看 


到 CME 的 扩展 。 


图 12 光 流 法 检测 CME 的 结果 ， 图 像 来 源 于 Colaninno “# (2006) P! 


ChinaXiv 合 作 期 刊 


Fig. 12 CMES detection results from an optical flow method, Image from Colaninno et al. (2006) [21] 

与 Colaninno 等 人 应 用 全 局 优化 的 光 流 算法 不 同 ，Gissot 等 人 ”提出 了 一 种 利用 梯度 估 
计 的 局 部 参数 化 光 流 方法 检测 CME， 通 过 Lucas-Kanade 算法 计算 相 邻 两 帧 图 像 在 某 一 时 间 
段 每 个 像素 点 位 置 的 移动 ， 即 速度 矢量 。 而 该 光 流 算法 有 一 定 的 约束 条 件 ， 仅 适应 于 小 的 目 
标 运 动 ， 当 目标 运动 速度 较 大 时 ， 即 帧 间 运 动 较 大 ， 算 法 误差 会 比较 大 ， 因 此 在 运动 分 析 中 
引入 基于 人 金字塔 的 多 分 辩 率 约束 ， 即 通过 二 次 采样 缩小 图 像 的 尺寸 , 同时 在 二 次 采样 前 对 图 
像 应 用 低 通 滤波 器 以 减 小 强度 变化 对 运动 估计 的 影响 , 提高 算法 的 鲁 棒 性 。 该 方法 通过 提出 
一 种 可 靠 运动 区 域 的 提取 策略 , 得 到 一 个 密集 的 速度 场 , 并 在 一 系列 连续 日 狗 图 像 上 进行 测 
试 来 检测 和 跟踪 相 邻 两 帧 图 像 之 间 的 运动 信息 。 

利用 光 流 法 检测 CME 运动 目标 的 基本 思想 在 于 给 图 像 中 的 每 个 像素 点 赋予 一 个 速度 矢 
量 ， 形 成 一 个 矢量 场 , 根据 各 个 像素 点 的 速度 矢量 特征 对 图 像 进行 动态 分 析 。 如 果 图 像 中 没 
有 运动 目标 ， 则 光 流 矢量 在 整个 图 像 区 域 是 连续 变化 的 ; 当 图 像 中 有 运动 物体 时 ， 目 标 和 背 
景 存在 着 相对 运动 ， 两 者 形成 的 速度 矢量 必然 不 一 样 ， 如 此 便 可 以 计算 出 图 像 中 CME 的 位 
置 。 
3.3.4 基于 学 习 的 方法 

传统 的 检测 方法 主要 通过 人 为 定义 特征 或 利用 设 定 简单 阐 值 进行 处 理 , 并 不 能 很 准确 地 
检测 出 CME 现象 ， 而 且 CME 通常 具有 明亮 的 结构 和 丰富 的 纹理 特征 ， 仅 使 用 单一 特征 识 
别 CME 区 域 也 不 够 准确 。 随 着 近年 来 机 器 学 习 与 深度 学 习 技 术 在 计算 机 视觉 与 图 像 处 理 中 
的 应 用 越 来 越 广泛 已， 机 器 学 习 中 有 很 多 效果 良好 的 分 类 器 被 用 来 解决 各 种 分 类 和 检测 
问题 。 同 时 CME 检测 是 一 个 二 分 类 问题 ， 即 识别 图 像 中 是 否 存在 CME， 因 此 基于 学 习 的 
方法 可 应 用 于 CME 检测 。 

由 于 CME 在 形状 和 空间 尺度 上 变化 很 大 ， 因 此 仪 使 用 单一 分 类 器 是 不 够 的 。 在 机 器 学 
习 领 域 , Yoav Freund 和 Schapire 提出 的 AdaBoost 算法 中 能 够 集成 同一 个 训练 集训 练 的 不 同 
弱 分 类 器 , 构成 一 个 强 分 类 器 , 该 强 分 类 器 的 性 能 超过 集成 的 任 一 弱 分 类 器 。 基 于 AdaBoost, 
Zhang 等 人 请 就 提出 了 一 种 CME 图 像 分 类 算法 ， 与 上 述 大 多 数 方法 在 极 坐标 变换 下 的 图 像 
中 检测 CME 不 同 ， 该 方法 直接 将 扇形 区 域 的 差分 图 像 分 割 成 块 ， 并 分 析 图 像 中 最 亮 块 的 灰 
度 。 通 过 实验 表明 : 切片 大 小 、 灰 度 阔 值 和 亮点 分 数 闵 值 都 会 影响 检测 结果 ， 因 此 该 方法 通 
过 获得 学 习 样 本 并 为 其 设计 分 类 模型 以 得 到 这 三 种 影响 因子 的 最 佳 参 数组 来 设计 弱 分 类 器 。 


= 最 后 使 用 AdaBoost 组 合 这 些 弱 分 类 器 ， 得 到 最 终 的 强 分 类 器 来 分 类 CME. 
O Hj Zhang 等 人 类 似 ，Yin 等 人 的 也 认为 CME 对 应 图 像 中 的 最 亮 块 ， 因 此 将 CME 的 检 


测 建 模 为 差分 图 像 中 最 亮 块 的 分 类 问题 ， 同 样 提 出 了 一 种 基于 AdaBoost 分 类 CME 图 像 的 
AE. 不同 的 是 该 算法 将 差分 图 像 进行 了 极 坐 标 变换 ,提取 变换 后 图 像 最 亮 块 的 灰 度 , 纹理 
All HOG 特征 ， 采 用 多 特征 融合 的 方式 构造 分 类 器 。 由 于 CME 具有 不 同 的 外 观 ， 灰 度 和 级 
理 特征 ， 因 此 选择 合适 特征 和 分 类 器 非常 重要 ， 而 决策 树 可 以 自动 实现 特征 选择 和 分 类 ， 因 
而 将 其 设计 为 弱 分 类 器 ， 同 时 AdaBoost 作为 集成 分 类 器 可 以 提高 单个 分 类 器 的 分 类 能 力 。 
于 集成 学 习 的 群体 智慧 和 决策 树 的 特征 选择 能 力 , 所 提出 的 基于 多 特征 的 检测 算法 的 集成 
以 获得 更 好 的 检测 效果 。 由 于 该 方法 主要 考虑 最 亮 块 区 域 , 而 忽略 了 上 暗 腔 区 域 及 其 他 类 似 
CME 的 太阳 结构 体 ， 容 易 漏 掉 比 较 微 弱 的 CME 事件 。 

针对 分 类 器 收敛 速度 慢 的 问题 , Zhang 等 人 ”还 提出 了 一 种 基于 极限 学 习 机 ELM) P?! 
的 新 检测 算法 ， 该 方法 通过 帧 间 差 分 、 极 坐标 变换 、 滤 波 去 噪 完成 预 处 理 。 与 上 述 他 们 提出 
的 方法 类 似 , 同样 基于 块 分 割 图 像 的 方法 , 不 同 的 是 之 前 提出 的 方法 是 将 最 亮 的 部 分 分 割 成 
块 ， 而 该 算法 是 利用 固定 大 小 的 块 遍历 预 处 理 后 的 整个 图 像 。 块 尺寸 的 选择 如 下 : 若 块 中 包 
含 的 所 有 像素 的 平均 灰 度 值 高 于 指定 闵 值 , 该 块 被 视 为 亮 块 , 通过 计算 每 个 块 大 小 下 的 分 类 
准确 率 (正确 检测 样本 与 所 有 样本 之 比 ) 选择 最 高 准确 率 的 块 尺寸 分 割 图 像 。 同 时 使 用 块 中 


=| 
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灰 度 和 纹理 特征 组 成 特征 向 量 并 构建 基于 ELM 的 分 类 器 。 该 方法 的 优点 是 基于 ELM 的 检 
测 方法 可 以 选择 合适 的 特征 ,收敛 速度 快 ， 同 时 通过 使 用 时 空 连续 性 的 决策 规则 ， 可 以 去 除 
大 多 数 具 有 类 似 CME 的 太阳 结构 体 ， 提 高 了 分 类 效果 。 

上 述 几 种 基于 机 器 学 习 的 方法 都 将 图 像 中 的 最 亮 块 视 为 CME， 因 此 都 是 对 最 亮 块 进行 
检测 和 分 类 ， 同 时 选择 的 分 类 器 也 各 有 特点 。 但 是 这 些 方法 都 是 通过 人 为 定义 特征 的 提取 ， 
HT CME 具有 多 种 特征 ， 人 工 选择 的 特征 不 一 定 能 得 到 好 的 检测 效果 。 而 深度 学 习 具 有 强 
大 的 学 习 能 力 和 特征 表达 能 力 ,会 自动 提取 特征 进行 分 类 。 姚 海 等 人 中 就 提出 了 一 种 基于 卷 
积 神经 网 络 (Convolutional Neural Networks，CNN) 的 CME 检测 方法 , CNN 采用 有 监督 的 学 
习 方 式 , 使 用 海量 具有 标签 的 训练 样本 ,结合 反 向 传播 算法 来 更 新 权重 ,训练 出 一 个 具有 分 
类 能 力 的 网 络 。 该 方法 自动 提取 适应 于 CME 检测 的 图 像 特 征 并 建立 检测 模型 ， 网 络 结构 采 
包含 卷 积 层 , 池 化 层 和 全 连接 层 的 六 层 卷 积 神经 网 络 。 为 了 网 络 模 型 能 够 在 训练 阶段 学 习 
到 较 好 的 参数 进行 分 类 任务 ， 训 练 采 用 明显 的 CME 图 像 数 据 集 ; 同时 以 此 参数 为 基准 在 微 
调 阶 段 采 用 微弱 的 CME 数据 集 进行 参数 调 优 ， 使 得 网 络 模型 能 获得 较 好 的 分 类 效果 。 相 较 
于 以 往 的 检测 方法 ， 基 于 卷 积 神经 网 络 的 分 类 方法 在 GPU 模式 下 处 理 图像 效 率 更 是 远 胜 于 
其 他 算法 ， 实 时 性 好 。 
目前 基于 学 习 方法 检测 CME 的 研究 还 不 太 多 ， 文 献 也 很 少 ， 与 上 述 传统 方法 相 比 ， 该 
类 方法 分 类 效果 好 ， 检 测 准确 率 高 ， 因 此 可 以 将 此 方法 作为 未 来 CME 检测 的 一 个 重要 研究 
方向 。 
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间 环 境 的 强烈 扰动 ， 对 通讯 、 导 航 、 航 天 器 等 高 技术 系统 可 能 造成 灾害 性 影响 ， 因 此 CME 
的 检测 对 于 预防 这 些 空间 灾害 性 空间 天 气 具有 非常 重要 的 作用 。 
目前 虽然 已 经 有 很 多 方法 应 用 到 CME 的 检测 上 ， 使 得 CME 的 检测 研究 取得 了 很 大 进 
展 ， 但 这 仍然 是 一 个 具有 挑战 性 的 问题 。 一 是 CME 具有 不 同 的 尺度 ， 大 小 和 形状 ， 这 使 得 
难以 选择 适当 的 特征 来 检测 所 有 CME; 二 是 CME 在 喷射 过 程 中 会 伴随 一 些 与 其 结构 非常 


c 相似 的 干扰 源 ， 因 此 很 难 将 CME 与 这 些 结构 区 分 开 来 ， 三 是 对 于 一 些小 的 ， 弱 的 CME, 
= 通常 会 发 生 误 检 与 漏 检 。 
© 而 随 着 近年 来 机 器 学 习 技 术 在 图 像 处 理 与 机 器 视觉 领域 获得 的 巨大 成 功 , 可 以 利用 机 器 


学 习 技术 替代 简单 的 阔 值 分 割 技 术 来 提高 CME 检测 效果 。 同 时 深度 神经 网 络 可 以 自动 地 提 
合适 的 特征 ,在 语音 识别 、 图 像 处 理 等 方面 都 取得 了 良好 的 效果 。 男 外 ， 现 阶段 大 数据 和 
高 性 能 计算 设备 都 已 完全 具备 ， 为 深度 学 习 在 CME 检测 中 的 应 用 提供 了 条 件 。 后 续 CME 

检测 的 主要 研究 方向 可 以 从 以 下 几 个 方面 开展 : 
C1) 继续 利用 机 器 学 习 的 方法 选择 合适 的 特征 以 及 分 类 器 ， 进 一 步 提高 CME 的 分 类 


效果 ; 


(2) 由 于 CME 的 爆发 过 程 是 一 个 动态 过 程 ， 因 此 可 以 将 基于 深度 学 习 的 运动 目标 检 
测 方 法 应 用 于 该 检测 ; 

(3) 可 以 认为 CME 的 检测 是 分 离 日 锡 图 像 中 前 景 运动 目标 的 过 程 ， 因 此 可 以 将 基于 
深度 学 习 的 前 景 检测 与 背景 减 除 方法 应 用 于 该 检测 。 

而 利用 基于 深度 学 习 方法 的 一 个 难点 和 重点 在 于 数据 集 的 制作 , 手工 标记 需要 大 量 的 人 
力 与 时 间 ， 同 时 对 一 些 弱 的 ， 小 的 CME 标记 时 可 能 会 有 误差 和 歧义 ， 因 此 还 需要 不 断 地 深 
入 研究 。 

总 而 言 之 , 随 着 人 工 智能 技术 的 发 展 和 进步 , 面向 大 数据 分 析 和 挖掘 的 各 类 基于 学 习 的 


技术 在 太阳 物理 学 研究 和 应 用 中 , 将 会 有 更 为 广泛 的 应 用 空间 和 应 用 前 景 。 因 此 基于 机 器 学 
习 和 深度 学 习 等 技术 可 以 很 好 地 应 用 到 CME 的 检测 上 ， 有 望 得 到 更 快速 ， 更 准确 的 检测 结 
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Research Progress On Detection Of Coronal 
Mass Ejection 


Guo Min, Li Jingtao, Shang Zhenhong, Liu Hui , Xian Xianggui , Yang Zhipeng 


( Kunming University of Science and Technology , School of Information Engineering and Automation , Kunming , 650500 , China ) 


Abstract: A Coronal Mass Ejection (CME) is a large-scale, severely explosive phenomenon 


that is the main solar outburst affecting the Earth. As this kind of explosion will cause serious 


interference to the terrestrial environment, detecting CMEs is of great significance for forecasting 


severe space weather. In order to more clearly sort out the existing CME detection methods, this 


paper analyzes and summarizes the typical methods. Firstly, CMEs and their characteristics are 


introduced. Then, CMEs detection methods are summarized and analyzed from two aspects: 


manual methods and automatic methods. Finally, problems of the current algorithms are discussed 


and the future research direction is proposed. 
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